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Abstract — This note presents an example of bilinear conser- 
vative system in an infinite dimensional Hilbert space for which 
approximate controllability in the Hilbert unit sphere holds for 
arbitrary small times. This situation is in contrast with the 
finite dimensional case and is due to the unboundedness of the 
drift operator. 

I. INTRODUCTION 

A. Control of quantum systems 

The state of a quantum system evolving in a Riemannian 
manifold fl is described by its wave function, a point ip in 
L 2 (il, C). When the system is submitted to an electric field 
(e.g., a laser), the time evolution of the wave is given, under 
the dipolar approximation and neglecting decoherence, by 
the Schrodinger equation: 

= (-A + V(x))1>(x, t) + u(t)W(x)^(x, t) (1) 

at 

where A is the Laplace-Beltrami operator on Q., V and W are 
real potential accounting for the properties of the free system 
and the control field respectively, while the real function of 
the time u accounts for the intensity of the laser. 
It is standard to rewrite the dynamics as 



^ = (4 + !l(t)% 



(2) 



where ip belongs to a separable Hilbert space H and (A, B) 
satisfies Assumption Q] 

Assumption 1: A and B are linear operators with domain 
D(A) and D(B) such that 

1) A is essentially skew-adjoint (possibly unbounded) 
with domain D(A); 

2) there exists an Hilbert basis (fa)keN of H made of 
eigenvectors of A. For every k, A fa = \\kfa', 

3) for every k in N, fa belongs to the domain D(B) of 
B; 

4) for every u in R, A + uB is essentially skew-adjoint 
on D(A) n D(B); 

5) B is essentially skew-adjoint. 

From Assumption Q] for every u,t in R, e t ( A + uB ) j s a un j. 
tary operator. By concatenation, for every piecewise constant 
function u, one can define the solution t H> T"^q of (01 with 
initial condition ip(0) = ?Ao- With extra regularity hypotheses 
(for instance: B bounded), it is possible to define T u for 
controls u not necessarily piecewise constant. 



A natural question, with many practical implications, is 
to determine the set of wave functions T^o that can be 
reached from a given initial condition -00 at a given time 
T when the control law u varies in the set U of admissible 
controls (here, U is the set of piecewise constant functions). 
This set {Tj.^>q, u £ U} is called the reachable set from iJjq 
at time T and is denoted TZr(ipo). For T £ (0, +oo), the 
reachable set from ijj at time smaller than T is 72. < t('0o) = 
U t <T^iC0o) and K <+oo (ip ) = U T>O ^ T (-0 O ). At this time, 
no definitive description of the reachable sets is known, with 
the exception of 7Z < + ca (ipo) for a few simple examples ([1], 
[2]) where U — L 2 ([0, +oo), R) and f2 is a bounded interval 
of R. Notice that even in this 1 -D framework, the results are 
far from obvious. 

Instead of looking at the complicated structure of 1Zt (ipo), 
one may consider its topological closure for a given norm. 
Many results of approximate controllability have been given 
in the last few years. We refer to [3], [4], [5], [6], [2], [7] 
and references therein for a description of the known theo- 
retical results concerning the existence of controls steering 
a given source to a given target. As proved in [4], [8], [9], 
approximate controllability is a generic property for systems 
of the type of (0. 

B. Temporal diameter 

Let a couple (A, B) of linear operators associated with an 
Hilbert basis (fa)keTSi of H made of eigenvectors of A be 
given. If 7£<+oo(0i) i s everywhere dense in the Hilbert unit 
sphere Sh of H, we define the temporal diameter of H for 
the control system © as 



p = inf jr>o 



K 



<T( 



The aim of this note is to give a positive answer (Theorem 
[Til l to the following question: "Does it exists a couple (A, B) 
satisfying Assumption [TJ such that p = 0?" 

The importance of this question goes beyond its purely 
mathematical aspect. A positive answer is a strong justifica- 
tion of the fact that we can neglect decoherence in (O, since 
decoherence has little effect in small times. 

C. Content of the paper 

The first part of this paper (Section HJi is a short summary 
of the finite dimensional case. It includes finite dimensional 
estimates (Section III-Bl i that are instrumental in our study. 



The second part (Section HlB presents basic infinite dimen- 
sional material. The third part presents an example (Section 
HV-Ab with zero temporal diameter. The technical computa- 
tion of this diameter (Sections HV-Bl and llV-Q combines the 
results of Sections IH-BI and IIII-DI 

II. FINITE DIMENSIONAL CASE 

In this Section, we concentrate on the case where H is 
finite dimensional. This case has been extensively studied 
([10], [11]). In this case, a classical choice for the set of 
admissible controls is the set L / 1 oc ([0, +00), R) of locally 
integrable functions in R.Other choices could be to restrict 
to the set of piecewise constant functions or to extend to the 
set of Radon measures. By continuity of the input-output 
mapping, and since we are only concerned with the closure 
of the reachable sets, all this choice are equivalent. 

A. Bilinear systems in semi-simple Lie Groups 

The operators A and B can be seen as skew-hermitian 
matrices. Changing A into A — Trace(^4)/dim(i/) and 
B into B — Tracc(i3)/dim(i/) just induces a physically 
irrelevant phase shifts. 

Considering the resolvent, the original system in H = C N 
can be lift to SU(N), the group of unitary matrices of order 
N with determinant 1. Denoting with x the matrix of the 
linear operator tp H> Tftp, (ffjl turns into 

j t x= (A + uB)x = dR x (A + uB) (3) 

with initial condition x(0) = Ih where R is the right 
translation in SU(N): R x : y H> yx. 

The lift of (fJI from C N to SU(n) as the right invariant 
controlled system (01 allows us to use the rich structure 
of semi-simple compact Lie groups. In particular, the Lie 
algebra su(N) = Ti N SU (N) turns into an Euclidean space 
when endowed with the(opposite) of the bi-invariant negative 
definite Killing form K(x,y) = (N ± l)Tia,ce(xy). 

The major drawback is that we have now to consider the 
whole propagator x (that is, equation (f2]i with all possible 
initial conditions). 

A classical result states that the temporal diameter of 
SU(N) for the dynamic of (f3]i is finite for generic pairs 
(A,B). 

Proposition 1: If Lie(A, B) — su(7V), then there exists 
T > such that TZ^itp) = Sh for every ip in S#. In other 
words, p < +00. 

Proof: We refer for instance to [12], ■ 

For every admissible control u, define v : t H> L u(t)At 
and y : t>-> exp (~v(t)B) x. The dynamics of y is given by 

j t y = e-^ B Ae^ B y := dR Ad ^ (t)B A y (4) 

The adjoint mapping Ad is an isometry for the Killing norm. 
In other words, the derivative of y has a constant Killing 
norm (equal to the Killing norm of A). As a consequence, the 
temporal diameter is positive as soon as the torus {e KB \K € 
R does not fill SU(n). 

Proposition 2: When H is finite dimensional with dimen- 
sion larger than or equal to 2, p > 0. 



B. Some time estimates 

In this Section, we give an estimate of the time needed 
to steer an eigenstate of A to another. The proof relies on 
averaging techniques (see [13] for details). 

Definition!: Let H = C N and {A<- N \ #>) satisfy 
Assumption!]] A couple (j, k) is a non-degenerate transition 
of (iW.fiW) if i) (^,#)^) ^ and ii) X 3 - X k = 
A; — X m implies that {j. k} = {I, m} or to} n {j, k} = 
or (0 /7 SW0 m ) =0. 

Proposition 3: Let (A^ N \ B^ N ') be a couple of skew- 
Hermitian matrices of order N. Assume that A^ N > is di- 
agonal and that (1,2) is a non degenerate transition of 
(iW.fiW). We denote with X? N) the propagator of x' = 



(AW 



tSWw, Define T 
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and u* 



t 1 Y 



|A 2 - Ai| cos((A 2 - X\)t + 1 
every t 



- IA2-A1I 
Then, for every n in N, for 



\X 



u* /n 
(N) 



(t,0) 



/(C + 1)||BW|| 



< 



1 + 2.K1IBWI 



with I = /„ |M*(r)dr = 4, K = ± f* |u*(*)|dt, A/ f the 
skew-Hermitian matrix of order N which entries are all zero 
but the ones of index (1,2) and (2,1) equal to 7r6i2e 1{?i /4 
and 7r&2ie _1 */4 respectively and 



C = sup 



J T M*(r)e i (^- A '=) r d7 



sm 



Ail 



where A is the set of all pairs (j, k) in {1, . . . , A^} 2 such that 
b jk ^ and {j, k} n {1, 2} ^ and \Xj - X k \ £ Z|A 2 - X x \. 

Proof: This is a particular case (for u* : t n- | A2 — 
Ai| cos((A 2 -Xx)t + 4>) of the inequality (13) in [13]. ■ 

C. A technical computation 

In order to apply Proposition [3] we will have to find K 
and (j> such that e KM sends a given vector to another one. 

For a given (a, f3) T £ C 2 , we will need r and in R 
such that e M (a,(3) T is colinear to (0, 1) T with 

re^ 
re'^ 



M 

For every s in R, 

exp(fAf) 

I or every (a, /3) T £ C 2 , 
a 



cos(r<) e 1 ^ sin(rt) 
e 1 ^ sin(rf) cos(rt) 



acos(rt) + e 1< ^sin(rt)/3 



exp(tAf) 
n the fol 

2 + |,3| 2 = 1. There exist 9 in [0,tt/2], a x ,j3i in (-7r,7r] 



P I V (3cos(rt) — e 1 ^ , asin(sa) 
In the following, we assume without lost of generality that 



\ a . 
such that 



a = cos 6e lai 
p = sin(9e i/31 . 



(5) 
(6) 



2 

9 + rt 



I(tt). 



With these notations, acos(rt) + sin(rt) = if and 
only if 

' cos( ai ~^~ ) cos(fl-rt) = 

sin(" Q1 -^ jcos(fl + rt) = 

Since we are interested in small t, we will chose <\> and t 
such that 

[ C os f 2i=&=4) = o . . 
<^ V 2 / that is, 

[ cos(0 + ri) = 

III. INFINITE DIMENSIONAL TOOLS 

A. Notations 

Let (A, B) satisfy Assumption Q] For every N in N, we 
define Cn the linear space spanned by <pi, fa, ■ ■ . , 4>n an d 
-kn : H ^ H, the orthogonal projection onto Cn'- 

N 
k=l 

The compressions of order TV of v4 and B are the finite 
rank operators A^ N ' = knA^k and B( N > = ttnB\c n - 
The Galerkin approximation of (O at order N is the infinite 
dimensional system 

^x = A^x + u{t)B^x. (7) 

Since £jv is invariant by (jTj, one may also consider (0 as 
a finite-dimensional system, whose propagator is denoted by 

The operator \A\ is a positive self-adjoint operator. In the 
case where A is injective, we define for every k > the 
norm /c-norm: ||^||fe = |||A| fe -0||. 

B. Basic facts 

Proposition 4: Assume that H has infinite dimension. 
Then, for every in S H , {e KB ^ 0l K e R} 7^ S ff . 

Proof: By the spectral theorem, up to a unitary trans- 
formation, H — L 2 (fl, C) with il a set of cardinality larger 
than two and B is the multiplication by a purely imaginary 
function. Then, the set {e KB i(>o, K G R} is included in the 
set {tp E iJsuch that |V>0e)| = IV'o (a:)|for a. e.x € fl}. The 
latter set is not equal to Sh- ■ 

Proposition 5: If H has infinite dimension and A is 
bounded, then p > 0. 

Proof: For every u in U, define Y u : t i-> 
e -/ |«(-r)|dTSYu^ p Qr ever y m g ffi f or ever y j n N 5 

the mapping t H> Y t u ipo satisfies 

In particular, for every ijj Q , for every t, 



< 
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Consider now two points an d ^1 m such that 
the distance 6 between the two sets {e KB ip$,K G R} and 
{e KB ipi,K e R} is not zero (such a couple (tp ,tpi) exists 



from Proposition |4j. Then p > S/\\A\\, since ^1 ^ IZ^ipo 
for t<S/\\A\\. M 
Proposition 6: If has infinite dimension and B admits 
an eigenvector in the domain of A, then p > 0. 

Proof: For every w in U, define as above Y u : t H> 
e~ lo |u( T )l dT - B Tf". v ^ e an eigenvector of B associated 
with eigenvalue A. For every ipo in Sh, the mapping t M> 
Y t u ipo satisfies 

A( w ,y t > ) = (^,r t ^ ). 

Fix in S# such that |(«,^o)| 7^ K u >V'i)l- Then p > 

<5/||Aw|| since Vi K <t ipo for t < S/\\Av\\. ■ 
Remark 1: Propositions [4] and [6] are true also when H 
has finite dimension larger than or equal to two. 

C. The RAGE theorem 

We define C b the linear space spanned by the eigenvectors 
(if any) of B. We define H B = C B 

Theorem 7: Let N in N and t/jo in H B . Then there exists 
a sequence (fc n ) nS N with limit +00 such that ||7rjve fe " B ^o|| 
tends to zero as n tends to infinity. 

Proof: This is a weak-version of the celebrated RAGE- 
theorem, see [14, Theorem XI. 115]. ■ 

D. Weakly-coupled quantum systems 

Definition 2: Let k be a positive number and let (A, B) 
satisfy Assumption Q] and such that the spectrum of iA is 
purely discrete (Xk)k and tends to +00 Then (A,B) is k 
weakly-coupled if for every u\ <E R, D(\A + u\B\ k / 2 ) = 
D(\A\ k / 2 ) and there exists a constant m such that, for 
every V in D(\A\ k ), \^t(\A\ k ^B^)\ < c {AiB) \(\A\ k ^)\. 
The notion of weakly-coupled systems is closely related to 
the growth of the |A| fc ' 2 -norm (\A\ k ip,ip). For k = 1, this 
quantity is the expected value of the energy of the system. 

Proposition 8: Let (A, B) be /c-weakly-coupled. Then, 
for every ip G D(\A\ k / 2 ), K > 0, T > 0, and u in 
^([0,00)) for which \\u\\ L i < K, one has ||T£(^o)|| fc / 2 < 

e c ^ B ^||^o|U/2. " 
Proof: This is [15, Proposition 2]. ■ 
Proposition 9: Let fc and s be non-negative numbers with 
< s < k. Let (A, B) be k weakly-coupled Assume that 
there exists d > 0, < r < k such that \\Bip\\ < d||V'||r/2 
for every tp in D(\A\ r l 2 ). Then for every e > 0, K > 0, 
n G N, and (ipj)i<j< n in £>(|A| fe / 2 )" there exists 7Y G N 
such that for every piecewise constant function u 



\ U \\L^ 



< K 



|T?(^)-X(jv)(*.0>JV^||./ a 



< e, 



for every t > and j = 1, . . . , n. 

Proof: This is [15, Proposition 4]. ■ 

Remark 2: An interesting feature of Propositions [8] and [9] 
is the fact that the bound of the |A| fc / 2 norm of the solution of 
(f2J or the bound on the error between the infinite dimensional 
system and its finite dimensional approximation only depend 
on the L 1 norm of the control, not on the time. 

Proposition 10: Let [A, B) be fc-weakly coupled for some 
k > 0. Then, for every in Sh, for every T > 0, 

{e K \JfeR}c^j, 



Proof: Fix K in R and e > 0. For every tj > 0, 
consider the control u^, constant equal to K/rj on the time 
interval [0, i]] and equal to zero elsewhere. By Proposition [9] 
there exists N such that ||T""(^ ) - XfyJt,0)%ifipo\\ < e 
for every i > 0. The classical theory of ODE ensures that 
X^(r],0)irNijjo tends to e KB irN"ipo as 77 goes to zero. ■ 

IV. AN EXAMPLE OF APPROXIMATE 
CONTROLLABILITY IN TIME ZERO 

A. A toy model 

We consider the following bilinear control system 

= -|A| a V> + u(i;)cos0V (8) 

where a is a real constant, £1 = R/27T is the one dimensional 
torus, H = L 2 (Cl, C) and A is the Laplace-Beltrami operator 
on f2. 

Remark 3: A realistic (and widely used) model for a 
rotating molecule is (JHJ with a = 1. For a^l, the presented 
example is purely academic. 

The self-adjoint operator — A has purely discrete spectrum 
{k 2 ,k € N}. All its eigenvalues are double but zero which 
is simple. The eigenvalue zero is associated with the constant 
functions. The eigenvalue k 2 for k > is associated with the 
two eigenfunctions 9 i-> cos(fc#) and i-> ^= sin(fc#). 
The Hilbert space H = L 2 (Vl,C) splits in two subspaces 
H e and H Q , the spaces of even and odd functions of H 
respectively. The spaces H e and H Q are stable under the 
dynamics of ©, hence no global controllability is to be 
expected in H. 

We consider the restriction of <[8J to the space H a . The 
function <pt : 9 1— > sin(fc0)/y / 7r is an eigenvector of 
the skew-adjoint operator A — i|Ai#J* associated with 
eigenvalue \k 2a . The familly {4>k)k&i is an Hilbert basis 
of H a . Here, B is the restriction to H of the multiplication 
by — 1003(6*). The skew-adjoint operators B is bounded and 
has no eigenvalue: Cb = H . For every j, k, ((pj, B<pu) = 
if j = k or \j - k\ > 2, and (cp^Bcpj+i) = -i/2. 

Theorem 11: If a > 5/2, then for every ipo,ipi in the 
Hilbert unit sphere of H , for every e > 0, for every T > 0, 
there exists a piecewise constant function u : [0, T] — > R 
such that ||T^-0o — V'lll < e - In other words, for ((H), if 
a > 5/2 then p = 0. 

B. Some time estimates 

Proposition 12: Assume a > 5/2. Let No large enough, 
P > No and ipo, ipi be in Cp n £]v _i such that ||V>o|| = 
1 1 '01 1| = 1- Then, for every s > 0, there exists u such that 
(|H) steers to an e-neighbourhood of ip\ in time less than 

604 1 2tt 

a 2 e{2a - 5) (AT - i)2«-4 + j^p- 
Proof: From Proposition (9j there exists iVi in N such 
that, for every u, \\u\\ L i < 8P implies that \\X^ Ni ^(t, 0)tpj — 
TV'jll — ^ or 3 = The problem is now to give an 
upper bound on the time needed to steer tp to ipi with system 

^ = A ^)^ + u(t)B^^. 



The idea is to find a control u\ that steers ipo to a neighbor- 
hood of e 1 ®°(j>N for some i?q m R- m time To and, similarly, 
a control 7/2 that steers ip\ to a neighborhood of e ^jvo for 
some $1 in [^0^0 + 27r) in time Ti. The final control is 
the concatenation of control u\, of control during time 
|$i — $o|/Aat and finally of the time -reverse of ui (which 
steers e 101 (pN o to a neighborhood of ip{). The total duration 
is Tq + |i?x — $o|/Ajv What remains to do is to give an 
estimate for To (the same computation gives a similar result 
for T x ). 

Let us describe intuitively our method. To induce the 
transition from ip to (pN a , we will first put the L 2 mass 
of the P th level of ipo to level (P — 1), without changing 
the modulus of the other coordinates of ipo- Then, we put all 
the L 2 mass of the (P — l) th level to level P — 2, and so 
on until all the mass is concentrated, up to a small error, on 
level Nq. 

We give now a formal description of the above method. 
Write ipo = 2~2n=n x n4>n- For N = No, ■ ■ ■ , P, we define 

u* :tH (Aat+i - Aw) cos(|A w+1 - A w |i + <^v) 



9m = arctan 



FjvI 



ifjv = 4- 



2(9 



iV 



j=iV+l Fj 



Jo 



((t)e i(((JV+2)^-(iV+l) 2 «)i dt 



sin ^7r 



(Af+2) 2 °-(jV+l) 2 



TN 



N e 
AN 2 



,n N 



(n+i) 2 <*-n 2 
4(l + 2K N \\B^\\)(C N 



l)\\B( N > 



2 A 



N+l 



Ajvl 



< 



8tt(1 



2K N \\B\\)(C N 



D\\B\ 



We proceed by induction on TV from P to No to steer ipo = 

ip^to^- 1 such that IKI-ttp-iJV'o _1 ||+lkjv -iVf < 

such that 

2 

At 



?yp. If ipo is constructed, we steer ipf? to V'o 

IKi-^^IP + ll^o-iV^ll^i-fi-E 



stepP-iVo, is in an J^n 2i]jv < e/2 neighborhood 

of e I1?0 (/>Ar for some i?o in R- 

This construction is done, following Proposition [3] by us- 
ing the control t M> u*(t)/njv during time tat. The L 1 norm 
of the control (equal to Kn) is chosen (using the computation 
of Section Hl-Q in such a way that exp (K^M*) sends the 
mass of the (N + l) th eigenstate to the N th one, that is a 
rotation of angle On- Since, for every N, bw,N+i = — i/2, 
the matrix ilP appearing in Proposition's a block diagonal 
matrix with every block equal to zero, but the one in position 
N, N + 1 equal to 

-ie^ N 7r/<' 
e -i VJV7I y 8 o 

The phase <pN is equal to — + tt where a^r and 
are the phase of coordinates N+l and N of ipo +1 - 

Notice that, while Proposition [3] is stated for n an integer, 
tin is not an integer in general. Nevertheless we can apply 



Proposition [3] with the integer part [njvj in order to ensure 
that the transformation is done with an error less than r\ n- 
Changing the integer ["ivj to the real number % does not 
change the bound on the error for more than an factor two 
for large rijv. 

With our definition of A and B, X N = N 2a and ||B|| = 
\/2/2. Straightforward computations yield, for No large 
enough, 

C N + 1<—, K N <4: t n < 1,1 



an 

and, for Nq large enough, 

p 

To < 



a 2 N eN 2a ~ A ' 



< 



< 



< 



151 



E 



u N=N 

151 2 
a 2 7V e(2a - 5) (No - 1 
302 1 



,2a-5 



a 2 e(2a-5) (iV -l) 2Q - 4 ' 

Finally, for No large enough, the total time needed to steer 
ipo to an e-neighborhood of ipi is less than 

604 1 2tt 



N 2a ' 



a 2 e{2a - 5) {N - l) 2 "" 4 
which concludes the proof of Proposition [12] ■ 

C. Proof of Theorem [77] 

Proof: [of Theorem QT) Let ^o> V'l m &H> £ > 0. To 
prove TheoremQTJ it enough to prove that the system ([8]l can 
approximately steer one point of {e KB ipo, K G R} to an e- 
neighborhood of one point of {e KB t/ji, K £ R} in arbitrary 
small time. 

Since a > 5/2, N^ am, 
tends to infinity. Define No in N such that 

604 1 2tt 



4 ~ 2q and Nn 2a tend to zero as N 



1) 



N 2a 



< il- 



ex 2 e (2a - 5) (No 

By the RAGE-Theorem (Theorem [7j, there exist Kq and K\ 
in R such that || 7Tjv e^^Vo II < £ and ||7Tjv o e ifl - B '0i|| < e. 
There exists P in N such that ||(1 — 7rp)e i: ' :oS ?/'o|| < £ an d 
||(l-7r P )e^i|| <e. 

By Proposition[T2l there exists a control that steers np(l — 
^N -i)ipo to an e-neighborhood of np(l — 7^0-1)^1 in time 
less than 77. This concludes the proof of Theorem QTJ ■ 

V. CONCLUSIONS AND FUTURE WORKS 

A. Conclusions 

The question of the minimal time needed to steer a 
quantum system from an arbitrary source to (a neighborhood 
of) an arbitrary target is of great importance in practice. 
This note presents a simple example of bilinear conservative 
control system in an infinite dimensional Hilbert space for 
which approximate controllability in the Hilbert unit sphere 
holds for arbitrary small time. 



B. Future Works 

At this time, we have no simple criterion to decide in 
the general case (unbounded drift operator A and control 
operator B without eigenvalue) whether a system of type 
(fJJ has zero temporal diameter. New methods will likely be 
needed for further investigations. 
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